Estimating e'/e in the standard model 
independently of Im X t 
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The CP- violating parameter e'/e is estimated in a novel way by including the explicit compu- 
tation of e in the ratio as opposed to the usual procedure of taking its value from the experiments. 
This approach has the advantage of being independent from the determination of the CKM pa- 
rameters Im At and of showing more directly the dependence on the long-distance parameters. The 
matrix elements are taken from the chiral quark model approach, the parameters of which have been 
recently determined by an updated fit of the CP-conserving amplitudes entering the AI = 1/2 rule. 
By sampling the ranges of the experimental inputs according to a normal distribution and those of 
the theoretical uncertainties according to a flat one, it is found that 

e'/e= (2.5 ±0.9) x 10~ 3 , 

in agreement with the current experimental data. A more conservative estimate 

0.9 x 10~ 3 < e'/e < 5.8 x 10~ 3 

is found by considering the smallest and the largest values obtained by taking a flat distribution for 
all inputs. 



In the standard model, the CP-violating parameter e'/e can be in principle different from zero because the 3x3 
Cabibbo-Kobayashi-Maskawa (CKM) matrix Vy, which appears in the weak charged currents of the quark mass 
eigenstate, is in general complex. On the other hand, in other models like the superweak theory [Q, the only source 
of CP violation resides in the K°-K° mixing, and e' vanishes. Establishing the precise value of e' is therefore of great 
importance. 

Experimentally e'/e is extracted, by collecting Kt, and K§ decays into pairs of ir° and 7T , from the relation 



Re e'/e ~ |?7 + _/?7oo| 2 - 1 /6 , 



(1) 



and the determination of r/oo and which are, respectively, the ratio of the amplitudes for Kt, 

Kg — > 7r°7r° and Kt, — > 7r + 7r _ over K$ — > 7r + 7r~. 

With the announcement last year of the preliminary result from the KTeV collaboration (FNAL) Q based on data 
collected in 1996-97, and from the NA48 collaboration (CERN) J| based on data collected in 1997-98, the long- 
standing issue of whether e' vanishes or not seems to be settled. By computing the average (see Fig. 3 below) among 
the two 1992 experiments (NA31 [f§ and E731 H) and the preliminary data of KTeV and NA48 one obtains 



Re e/e' = (1.9 ±0.46) x 10" 



(2) 



The error in (g) has been inflated according to the Particle Data Group procedure for combining results with sub- 
stantially different central values. The value and range in can be considered the current experimental result. Such 
a result will be further improved by the complete run and full data analysis from KTeV and NA48 and the first data 
from KLOE (Frascati) thus achieving an uncertainty of only few parts in 10 4 . 

From the theoretical point of view, e' , which parameterizes direct CP violation in the decays (for a review see, e.g., 
is computed as 
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where ImAt = lmVtdV t *, 



1 



Ho = — V y> Re (Qi) (1 - fW) , (4) 

COS On — 

n 2 = ^E^ Re ^)^ ( 5 ) 

i 

and (<3i)i = (27r, Z|<5i|-?f). The phases come from final state interactions. The explicit forms of AS = 1 four-quark 
operators Q% can be found, for instance, in The term 0,,+,;' is the isospin breaking (for m u ^ m^) contribution of 
the mixing of ir with r\ and 77'. In (||), the phase <\> = tt/2 + Sq — 82 — 9 e = (0 ± 4)° has been taken as vanishing , and 
CPT is assumed to hold. The Wilson coefficients yt are known to the next-to-leading (NLO) order in a s and a w JlO| . 
Notice the explicit final-state-interaction phases Si which comes from writing in (||) and (0) the absolute values of the 
amplitudes in term of their dispersive parts. 

The parameter e, which parameterizes indirect CP violation, is given by 

e = lm\ t C e B K {ReA c [t7iSo(x c ) - rj 3 Si(x c ,x t )} - Re A^So^t)} (6) 

where A c = V c dV* s , Si(x) are the Inami-Lim functions %i — m i/ m w anc ^ G £ ^ s a constant equal to 

G%f^mKm^ v /(3V2TT 2 ArriLs)- Br is the bag parameter to be determined in the long-distance estimate of the 
hadronic matrix element for the AS* = 2 transition. In (^), a term of higher order in A = V us is neglected. The 
parameters rji , which encode the renormalization group running, are known to the NLO order . 

In the usual approach — followed by all current estimates — the ratio e'/e is computed by dividing (Q) by the ex- 
perimental value of e. The allowed values for the CKM combination Im A t in front of (|^) are then either taken from 
independent analysis of the unitarity triangle |l3|-|T5|] or consistently computed by means of a determination of e in 
the same model (which includes the long-distance matrix element of the AS" = 2 amplitudes) [p^|-|l8| . 

This two-step procedure can be by-passed by computing in a given model e' and e and by taking directly their 
ratio. Because both parameters are proportional to ImAt, this quantity, which is an important source of uncertainty, 



simplifies in the ratio. This approach has also been advocated in 19 16|] without, however, any numerical analysis. 
Of course, one trades the dependence on the Wolfcnstcin parameter 77 for that on p in the factor Re At = — A(l — 
A 2 /2)Vj ) (l — p) that enters in (^J). However, this dependence is weaker and p = (1 — A 2 /2)p can be determined 
reasonably well in the unitarity triangle independently of kaon physics. Another advantage of this procedure is that 
the dependence of the final result on many of the input variables is made more transparent because the computation 
is not separated into the determination of Im At on the one hand and e'/e on the other. However, as we shall see, the 
uncertainty in the final estimate is about the same size as that obtained by the usual method, and it is still dominated 
by the uncertainty in the hadronic matrix elements. 

In this work I use this second approach and apply it to the estimate of e'/e in the standard model. 

The hadronic matrix elements for all the relevant quark operators Qi_io and the parameter Bk are computed in 



the chiral quark model (xQM) pO| , pT[ at 0(p 4 ) in the chiral expansion along the lines presented in |18|. I use this 
approach because it is the only theoretical model in good agreement with the experimental data. 
As we shall see, the final result is consistent with that of the xQM obtained by the usual method. 
In the xQM approach, there are three model-dependent parameters M, (qq) = ((qq)) 1 / 3 and (a s GG/ir) which are 
fixed by means of a fit of the CP-conserving amplitudes in the AI =1/2 selection rule of K — > irir decays. This 
fit has been re-done recently in [ p2[ by updating those short-distance inputs that have become better known in the 
meantime. The matching of the hadronic matrix elements to the Wilson coefficients is done at the scale [i = 0.8 
GeV. The amplitudes A(K° -2n,I = 0) and A(K° - 2tt, I = 2) computed in the model are compared with their 
experimental values by allowing at most a ± 20% uncertainty. The values for M, (qq) and (a s GG/n) found are 
those in Table I. The values for the two condensates vary according to the 75-scheme of dimensional regularization, 
't Hooft-Veltman (HV) or Naive Dimensional Regularization (NDR), used in the computation of the NLO Wilson 
coefficients and hadronic matrix elements. 

Having thus determined the model-dependent parameters, it is possible to compute e'/e. 

Notice that the uncertainty in £l v + n > ]23| ] affects the final estimate only marginally since any change in this input 
implies a change in the fit to the AI = 1/2 rule with the net effect of recovering the original prediction. In fact, any 
change in 0^+^' is anti-correlated to the value of (a s GG/n) oobtained in the fit, the variation of which compensates 
in e'/e the original change of £!,,+,,<. Changing the value of m c affects the value of e'/e below the 10% level. 

The value of e'/e depends — besides other given parameters like meson masses and decay constants, and the matching 
scale p, that we keep fixed — on eight input parameters: M, (qq), (a s GG/ir), m s , Aqcd, mt, Vcb and p. The range for 
p is that obtained in the unitarity triangle without the bounds from the experimental value of e. See Table I for a 
summary of the ranges used. For central values of all input parameters, e'/e is equal to 2.5 x 10~ 3 in both 75-schemes. 
This value corresponds to the determination Bk = 1-0. By comparison with the experimental value of e, it is then 
found that Im A ( = 1.3 x 10~ 4 and 1.2 x 10 -4 , respectively, in the HV and NDR scheme. 
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parameter 

yi QCD 
[V<*| 
P 

m t (m t ) 

T 7+ T 7 / 

m a (1 GeV) 

m 

(asGG/n) 
M 



value 

340 ± 40 MeV 

0.0405 ± 0.0015 

242+ 070 
U.Z4Z_ 045 

165 ± 5 GeV 
0.25 ±0.10 
150 ± 25 MeV 
HV: (-210 <-> -260 MeV) 3 , NDR: (-230 
HV: (325 335 MeV) 4 , NDR: (327 <- 
HV: 180 ^ 220 MeV, NDR: 185 «-> 



<-> -260 MeV) 3 
■ 340 MeV) 4 
210 MeV 



TABLE I. Table of the numerical values of the most relevant input parameters used in the present analysis. 

In order to estimate the error, a set of values for e'/e is generated by varying the parameters M, (a s GG /it) , (qq) 
in a flat scanning of their given ranges, while the others according to a normal distribution. The difference in the 
values of e'/e found in the HV and NDR schemes is rather small in the xQM approach and it can be considered as 
part of the overall theoretical uncertainty by joining the values found in the two schemes. The set of values of e'/e 
thus found is amenable of statistical analysis. For a given set, a distribution is obtained by collecting the values of 
e'/e in bins of a given range. This is shown in Fig. |lj. The final distribution is partially skewed, with more values 
closer to the lower end but a longer tail toward larger values. 



0.7 1.2 1.7 2.3 2.8 3.3 3.9 4.2 4.4 4.9 5.4 




0.7 1.2 1.7 2.3 2.8 3.3 3.9 4.2 4.4 4.9 5.4 



FIG. 1. Distribution of values of e'/e . Normalized bins versus the central values (xlO 3 ) of each bin. 

However, because the skewness in the distribution is less than one (and the standard deviation of the variance much 
smaller than the variance itself), the mean and the standard deviation are a good estimate of the central value and 
the dispersion of values around it. The statistical analysis therefore yields 

e'/e = (2.5 ±0.9) x 1(T 3 , (7) 

which is the best estimate quoted in the abstract. 

The lcr uncertainty of the estimate (0) should be considered together with the range given by the largest and the 
smallest values of e'/e when all input parameters are varied in a flat scanning within their errors. In this case, it is 
found that 

0.9 x 1(T 3 < e'/e < 5.8 x 1CT 3 , (8) 

which is the more conservative result. 

It is also possible to gain some insight on the dependence of these results on the input parameters. For instance, in 
Fig. ||, the generated values of e'/e are plotted against, respectively, (qq) and m s while randomly varying the other 
parameters. In particular, larger values of (qq) correspond to larger values of e'/e. 

There is no strong dependence on m s or (a s GG /tt) because m s only enters in the NLO corrections and the gluon 
consensate does not contribute to the leading penguin operators. However, these two parameters play a role in the 
determination of e where the final values of -Br- depends on m s and (a s GG /it) . Given (||), e'/e is inversely proportional 
to the value of Bk ■ 
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FIG. 2. Variations of e'/e as the quark condensate (qq) and m s are varied. Other inputs are randomly generated within 
their bounds. 



The results in (0) and (g) are depicted in Fig. 3 as yellow/light-gray and gray bands together with the most recent 
experimental determinations and the average in (2). As it can be seen from this figure, the result is in good agreement 
with the experimental data. 
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FIG. 3. Experiments vs. the present estimate of e'/e. The dashed line is the central value, the gray area spans the conservative 
range (^) and the yellow/light-gray area spans the la range in ((?]). 

As discussed in more details somewhere else Q, the rather large values of e'/e found by means of the xQM 
hadronic matrix elements, and the corresponding better agreement with the experiments, comes about because the 
fit of the Al = 1/2 rule — on which the approach is based — enhances the contribution of the isospin 1 = channel 
which, in turns, drives e'/e toward larger values. 



I thank my collaborators S. Bertolini and J. Eeg for discussions and F. Parodi and A. Stocchi for their help in the 
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